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ABSTRACT 

Long Gamma-ray Bursts (GRBs) signal the death of massive stars. The afterglow 
emission can be used to probe the progenitor ambient through a detailed study of 
the absorption pattern imprinted by the circumburst material as well as the host 
galaxy interstellar medium on the continuum spectrum. This has been done at optical 
wavelengths with impressive results. Similar studies can in principle be carried out 
in the X-ray band, allowing us to shed light on the material metallicity, composition 
and distance of the absorber. We start exploiting this route through high resolution 
spectroscopy XMM-Newton observations of GRB 090618. We find a high metallicity 
absorbing medium (Z £ 0.2 Zq) with possible enhancements of S and Ne with respect 
to the other elements (improving the fit at a level of > 3.4 a). Including the metallicity 
effects on the X-ray column density determination, the X-ray and optical evaluations 
of the absorption are in agreement for a Small Magellanic Cloud extinction curve. 

Key words: gamma-rays: bursts - X-rays: general - X-ray: ISM 



1 INTRODUCTION 

Tracking the evolution of metals from star-size scale to large- 
scale structures is a major step in understanding the evolu- 
tion of the Universe. Metals are produced in stellar interiors 
and then ejected in their environs through supernova (SN) 
explosions and stellar winds, thus enriching the interstellar 
medium (ISM) of their galaxy. So far, the study of metal en- 
richment in galaxies at high redshifts has been carried out 
mostly by observing galaxies along the line of sight of bright 
quasars. This technique is plagued by selection effects, e.g. 
radiation from quasars probes more probably halos of the in- 
tervening galaxies. Gamma-ray bursts (GRBs) are opening 
a completely new window in understanding the history of 
metals and galaxy formation, in particular at high redshift. 
This is because, for geometrical reasons, quasar sight lines 
should preferentially intersect the outer regions of the ISM 
in high— z galaxies. In contrast, GRBs originate within the 
densest regions of their host galaxies where massive stars 
are produced, probing denser environments (Prochaska et 
al. 2007). It is now recognized that long-duration GRBs are 
linked to collapse of massive stars, based on the association 
between (low— z) GRBs and (type Ic) core-collapse super- 



novae (Woosley & Bloom 2006 and references therein). In- 
deed, at least ~ 75% of GRBs with known redshift show 
intrinsic (i.e. at the redshift of the host galaxy) X-ray ab- 
sorption larger than 3 x 10 21 cm -2 (Campana et al. 2010; see 
also Stratta et al. 2005; Campana et al. 2006). In particular, 
more than 80% of GRBs show an intrinsic column density 
larger than 5 x 10 21 cm " 2 (Campana et al. 2010). Fur- 
thermore, the presence of absorption indicates a significant 
metal enrichment, because X-rays are effectively photoelec- 
trically absorbed only by metals. Metallicities larger than 
those derived from Damped Lyman Alpha (DLA) studies at 
the same redshift are also measured from optical studies of 
GRB. They indicate that the metal abundance can be as 
high as 10% of the solar one up to z — 6 (Savaglio 2006; 
Kawai et al 2006; Campana et al. 2007). Signatures for the 
presence of this material can be expected in the optical light 
curve of GRB afterglows showing up as flux enhancements 
in their light curves or as (variable) fine-structure transition 
lines in their spectra (Vreeswijk et al. 2007; D'Elia et al. 
2009a, 2009b) as well as Lya (Thone et al. 2011). The ob- 
servations of these features can set important constraints on 
the density and distance of the absorbing material located 
either in the star-forming region within which the progenitor 
formed or in the circumstellar environment of the progeni- 
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Table 1. XMM-Newton observation log. 



Instrument Start time Exp. time Net exp. time 





(hr after burst) 


(b) 


(s) 


MOS1 


13:48:09 (5.33) 


22615 




MOS2 


13:48:08 (5.33) 


22620 




pn 


14:10:31 (5.70) 


21034 


12695 


RGS1 


13:47:24 (5.32) 


22842 


19675 


RGS2 


13:47:32 (5.32) 


22789 


19651 



tor itself (Prochaska et al. 2006; Dessauges-Zavadsky et al. 
2006; D'Elia et al. 2009a, 2009b). 

Metal abundance derived from optical spectroscopy can 
be underestimated because some of the metals can be locked 
in dust grains, the composition of which can hardly be de- 
rived through optical spectroscopy. This is not the case of 
X-ray measurements: X-rays are photoelectrically absorbed 
and scattered by neutral and ionised atoms as well as by 
dust grains, providing a full account of the metal content. 
Indeed, the soft X-ray band is the best suited for metallicity 
diagnostics: absorption comes from the innermost shells (K) 
of elements from carbon (0.28 keV) to zinc (9.66 keV), with 
oxygen (0.54) and iron (7.11 keV) being the most prominent. 
Iron L-shell edges are also detectable. In the X-ray domain, 
however, relatively little progress has been achieved in the 
study of GRBs, mostly due to low statistics, to low spectral 
resolution and to the fact that some absorbing features are 
redshifted outside the instrument spectral band. A transient 
feature consistent with being a Fe absorption has been re- 
ported for GRB990705 and GRB011211 (Amati et al. 2000; 
Frontera et al. 2004) in the prompt emission phase. More re- 
cently, a decreasing total absorbing column density has been 
reported in one of the farthest GRB050904 (at z = 6.3) and 
this has been modeled, together with optical data, to unveil 
the high metallicity of the circumburst medium (Campana 
et al. 2007). None of these studies however have been able 
to shed light on the chemical properties of the progenitor 
because the relatively poor photon statistics did not permit 
a more detailed analysis of the data. 

Here we report on an XMM-Newton observation of the 
bright GRB 090618 started 5.3 hr after the GRB onset. In 
Section 2 we describe the GRB 090618 properties. In Section 
3 we describe XMM-Newton data extraction and analysis 
and in Section 4 we discuss our results. 



2 GRB 090618 

The Swift Burst Alert Telescope (BAT) detected and located 
GRB 090618 on June 18 2009, 08:28:29 UT (Schady et al. 
2009). The burst was bright with a 1 second 15-150 keV 
peak flux of 38.9 ± 0.8 ph cm~ 2 s _1 . The burst duration 
(as measured by the time including 90% of the total flux, 
T 90 ) in the 15-350 keV energy band is T 90 = 113.2 ± 0.6 s. 
The GRB monitor onboard Fermi also observed this burst 
finding a peak energy of E p = 155.5l}J'g keV and a 8-1000 
keV fluence of (2.70 ±0.06) x 10~ 4 erg cm 2 (McBreen 2009). 

The X-ray telescope (XRT) started observing GRB 
090618 124 s after the onset. The burst was initially very 
bright also in X-rays, with a peak rate of ~ 8,000 c s~\ 
and then decayed rapidly with a power-law slope of ~ t~ 6 . 
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Figure 2. Contour plot of metallicity of the absorbing medium 
intrinsic to the host galaxy (with fixed solar abundance pattern) 
vs. the medium column density. The cross indicates the best fit. 
Lines mark 1, 2 and 3 a confidence levels. 



The afterglow decay showed a break around 310 s flattening 
to t~ ' 7 . The light curve than breaks at txi = 5, 500 s to a 
steeper decay index t~ 1A , and then again at tx2 = 3.1 x 10 5 
s to a decay as i -1,9 (Schady, Baumgartner & Beardmore 
2009). 

The optical afterglow was detected by the P60 telescope 
approximately 2 minutes after the burst. The initial magni- 
tude was r' ~ 13.8 (Cenko 2009). The afterglow showed a 
decay and a rebrightening reaching a peak at 120 s with a 
first break around ~ 600 s (Li et al. 2009). The afterglow was 
detected in all the Swifb-XJVOT filters (Schady 2009) and 
the light curve evolution is different from what observed in 
the X-ray band (Schady, Baumgartner & Beardmore 2009). 
Malesani (2009) noted that at the position of the afterglow 
there is a bright host galaxy with r = 22.7 and i = 22.2. 
Cano et al. (2011) collected optical information on this burst 
and provided evidence for the presence of a supernova at late 
times. 

Given the afterglow brightness a spectrum has been ob- 
tained with different facilities. The 3-m Shane telescope at 
Lick observatory observed GRB 090618 about 20 minutes af- 
ter the burst. A number of strong absorption features iden- 
tified as Mg II, Mg I, and Fe II at a common redshift of 
z = 0.54 were observed (Cenko et al. 2009). Fatkhullin et al. 
(2009) confirmed the redshift based on the detection of an 
[Oil] line. 

At z = 0.54 and using standard cosmology the isotropic 
equivalent energy in the 8-1000 keV band is E iso = 2.0 x 10 53 
erg, making GRB 090618 one of the brightest (in terms of 
total energy emitted) burst. In addition, at z = 0.54 GRB 
090618 is one of the very few GRBs nearby and intrinsically 
very bright. 
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Figure 1. XMM-Newton spectrum of GRB 090618 and X-ray continuum model. The upper datasct (black) refers to EPIC/pn data. The 
middle spectra are RGS1 (red) and RGS2 (green) first order spectra. The lower spectra are RGS1 (blue) and RGS2 (magenta) second 
order spectra. The best fit power law model with intrinsic absorption with solar metallicity is shown with a solid line. In the lower panel 
the contribution to the x 2 of each channel (retaining the sign) is shown. 



Table 2. Spectral results. 



Model 


N H (z) 
(10 21 ) 


Metallicity 
(solar units) 


Power law 
photon index 


X 2 (dof) 


Fix. Metall. 
Free Metall. 
Free. Metall+Ne+S+Fe 


n o+O.l 

z -°-o.i 

5-tds 

12±S 


1 (fix.) 
> 0.2 
< 0.6 


I Q1+0.02 
±,a -0.02 

1 89+ - 02 

L - a -0.02 


1288.6 (1262) 
1286.8 (1261) 
1263.1 (1258) 



3 XMM-NEWTON OBSERVATIONS 
3.1 Data preparation 

XMM-Newton started observing GRB 090618 on June 18, 
2009 13:48:09.0 UT, 5.3 hr after the burst onset, during 
j-1.36 jjgpay p ar t of the afterglow light curve. The expo- 
sure time was 22.6 ks for the MOS1 and MOS2 instruments, 
21.0 ks for the pn instrument (all in full window mode), and 
22.8 ks for the RGS1 and RGS2 instruments (see Table 1). 
Given the source strength MOS data are piled-up and were 
not considered in the following analysis. Data were processed 



using SAS version 9.0.0, using the most updated calibration 
files available. Standard data screening criteria were applied 
in the extraction of scientific products. Data of the pn in- 
strument were filtered to a 0.3-10 keV rate over the entire 
instrument less than 100 c s -1 , resulting in 12.7 ks net ex- 
posure time. This filtering threshold has been selected as 
a trade-off between low background and exposure time. We 
did not filtered out RGS data. Despite its strength, the GRB 
afterglow emission pn data were not piled-up (less than 1% 
at the beginning of the observation when the GRB afterglow 
emission was stronger). RGS data, thanks to the dispersion 
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of the counts over a larger detector area, are not piled-up. 
The pn data were extracted from a large circular region with 
radius 1,600 pixel (80 arcsec, pattern to 4). We extracted 
85,000 counts for a mean count rate of 5.8 c s _1 . The back- 
ground was extracted on the same CCD of the source from a 
1,030 pixel radius circular region (in order to avoid the GRB 
streak) at ~ 5' from the source. The background accounted 
for less than 10% of the source counts. The response ma- 
trix and the ancillary response file were generated with the 
dedicated software tasks. 

RGS data were reduced using the rgsproc task. First 
and second order background subtracted spectra for the two 
instruments were retained and the appropriate response ma- 
trices were generated. RGS data were slightly filtered in or- 
der to avoid the strongest part of the solar flares, resulting 
in 19.5 ks. In total we collected 4125 (1579) for the first 
(second) order spectrum with RGS1 and 4381 (1504) for 
the first (second) order spectrum with RGS2, respectively. 
The background for these spectra is around 25 — 40% of the 
source flux. 

3.2 Data analysis 

The spectra are fit with a power law spectrum absorbed by 
a Galactic component and by an intrinsic (z = 0.54) com- 
ponent. In order to assess the contribution of the Galac- 
tic absorption we first considered HI maps which provide a 
weighted mean column density at the position of the GRB 
of 5.8 x 10 20 cm' 2 (Kalberla et al. 2005) and 6.5 x 10 20 cm~ 2 
(Dickey & Lockman 1990). To fit the absorption component 
we adopted the Wilms et al. (2000) model (tbabs within 
the X-ray fitting program XSPEC) and abundance pat- 
tern. To check the Galactic absorption we fit the data of 
a closcby (8' from the GRB) AGN (R.A.: 19h35m57.5 Dec: 
+78d21m27.0) with a power law model. The resulting col- 
umn density is 9.5±3;f x 10 20 cm~ 2 (90% confidence level). 
This is in line with the predictions by HI maps even if some- 
what on the high side. Despite its small influence, we fix the 
Galactic contribution to 6.5 x 10 20 cm~ 2 . We perform the 
first fit grouping the RGS data to 20 counts per spectral 
bin and the pn data to 50 counts per bin in order to have a 
good tracing of the continuum at high energies. A fit with 
the metallicity fixed to the solar values gives an intrinsic 
column density at z = 0.54 of N H (z) = 2.3^'J x 10 21 cm~ 2 
with a x 2 = 1288.6 with 1262 degrees of freedom (dof, i.e. 
a reduced x?ed = 1-02) and 29% null hypothesis probability 
(nhp). The power law photon index is Y = l^lt^g 2 , (see 
Table 2). We also checked that the column density did not 
decrease by splitting into two parts the entire observation 
with almost the same number of counts. 

Given the relatively large number of counts collected we 
leave free to vary the absorbing medium metallicity (keeping 
the solar abundance pattern). The fit do not improve signif- 
icantly with x 2 = 1286.8 with 1261 dof (xLj = 1-02) and 
30% nhp. An error search shows that the metallicity of the 
absorbing medium is ^ 0.2 Zq (see Fig. 2). The intrinsic col- 
umn density is N H (z) = 5-0±l.l x 10 21 cm~ 2 . For the lower 
limiting metallicity the column density is < 7.7 x 10 21 cm~ 2 . 

The number of counts is not high enough to afford an 
element by element analysis. To circumvent this we leave free 
the abundance of the single elements, one at a time, keeping 
the others all tied together. For many of them we just find 



Table 3. Abundance limits*. 



Element 


Abundance 


Element 


Abundance 




(solar units) 




(solar units) 


C 


< 4.9 


N 


< 4.9 


o 


< 0.6 


Ne 


0.4-0.7 


Mg 


< 0.7 


Si 


< 2.2 


S 


2.6 - 25 


Fe 


0.2 - 0.8 



* Upper limits are 90% confidence level for one parameter of 
interest (Ax 2 = 2.71). Intervals are 90% confidence level for one 
parameter of interest. These values were obtained by leaving free 
the abundances of each element, one by one keeping the other 
abundances tied together and free to vary. 



upper limits (see Table 3). In particular, a few elements show 
an abundance not consistent with zero, namely Ne, S and 
Fe. 

Motivated by this analysis we fit the same spectra leav- 
ing free the abundances of Ne, S and Fe and linking together 
all the other elements (mainly driven by Oxygen). The fit 
improves with a x 2 = 1263.1 with 1258 dof (x 2 cd = 1-00) 
and 45% nhp. Formally, an F-test gives a random proba- 
bility for the improvement of 3.3 x 10~ 5 , equivalent to a 
significance of 4.1 a. The new fit provides a column den- 
sity of Nh = l-2lo'g x 10 22 cm" 2 and just an upper limit 
on the metallicity of the remaining elements of < 0.6 Zq. 
The abundances of the three elements left free to vary to- 
gether lie within Ne = 0.4 - 1.8 Ne , S = 1.5 - 19 S Q , and 
Fe < 0.3 Feo, respectively. A contour plot of S and Ne abun- 
dances is shown in Fig. 3. We also searched for structures 
due to absorption of different ionization stages and/or to 
dust around the edges of the above elements at the known 
GRB redshift of z = 0.54, finding none. Given the presence 
of possible uncertainties in the Galactic column density de- 
termination, in the adopted solar abundance pattern we left 
free to vary the Galactic column density in a ±50% range 
around the selected value. We then repeat the fit with a 
constant metallicity and the fit with Ne, S and Fe free to 
vary. The fit with constant metallicity is characterized by 
a large Galactic column density (at the edge of the allowed 
distribution) and with a xLi = 1.02 (1260 dof). The fit with 
variable Ne, S and Fe abundances gives a Xrcd = 1-00 (1257 
dof). In this case the F-test gives a a random probability for 
the improvement of 6.0 x 10~ 4 , equivalent to a significance 
of 3.4cr. 

We also search for possible emission lines. We first 
search for an iron line looking for an unresolved emission line 
in the 6.4-6.9 keV (rest frame) energy interval. This range 
is covered by pn data only. We do not find any significant 
line in this range and we are able to put an upper limit (90% 
confidence) of 37 eV to any line. We also inspect the RGS 
data. One of the strongest (narrow) emission lines is OVIII 
at 654 eV (see, e.g., Bertone et al. 2010 for the list of line 
searched). We find a hint for the presence of this line at an 
energy of 432 eV (665 eV rest frame) even if its significance 
is between 90% and 99% (see Fig. 4). The line equivalent 
width is 4.8 eV. The velocity difference of the line energy 
with respect to the quoted redshift is on the blue side of 
the spectrum, possibly indicating an outflow, but amounts 
to 5, 000 ± 650 km s" 1 , which seems somewhat high. 
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Figure 3. Contour plot of the Ne vs. S abundances in solar units. 
The cross indicates the best fit. Lines mark f, 2 and 3 cr confi- 
dence levels. The spiky features present in the contour plot are 
not real and are likely due to numerical problems induced by 
the constrained interval within which the Galactic column den- 
sity lies. The real contour level should be obtained by joining the 
outermost values. 
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Figure 4. Contour plot (line energy vs. line count rate normal- 
ization) of the putative O VfH detected in the XMM-Newton 
data. Energies are in the observer frame. Lines mark 1, 2 and 3 
a confidence levels. 



4 DISCUSSION AND CONCLUSIONS 



ponent and an intrinsic component at the GRB redshift 
(z = 0.543). 

Assuming the Wilms et al. (2000) solar abundance 
patter, we estimate the mean metallicity of the absorbing 
medium close to the GRB. The metallicity of the medium 
is strongly covariant with the total amount of the absorb- 
ing column (see Fig. 2) and we are able to set a lower limit 
of 0.2 Zq (90% confidence level). The collapsar scenario 
for GRBs requires a massive helium star with rapidly rotat- 
ing core (Woosley 1993; MacFadyen & Woosley 1999). Low 
metallicities are requested by models in order to ensure that 
the rotationally induced mixing process produces a quasi- 
chemically homogeneous stellar evolution avoiding the spin- 
down of the stellar core (Izzard, Ramirez- Riuz & Tout 2004; 
Yoon & Langer 2005; Woosley & Heger 2006). This metal- 
licity threshold is expected to be around 0.1 — 0.3 Zq. If 
the medium probed by our observations has not been heav- 
ily enriched by the progenitor and is therefore indicative of 
the progenitor's metallicity, our result is (barely) consistent 
with model predictions. Given however the low redshift of 
the GRB and the brightness of the host galaxy it can also be 
expected that the mean metallicity of the medium is high. 

Analyzing the pn and RGS data we found that the ab- 
sorbing medium might be overabundant in Ne and S (see 
Fig. 3 and Table 3). Both best fit column densities are 
~ 6 x 10 17 cm" 2 . Clearly we do not have an idea of the 
distance of the absorbing medium nor of its geometry. Con- 
sidering spherical symmetry and a distance scale of 0.1 pc 
we have 0.01 Mq and 0.02 Mq mass of Ne and S, respec- 
tively, as well as 11 Mq of hydrogen. This sets the order 
of magnitude of the distance in the case of a medium en- 
riched by the progenitor. Larger distances naturally involve 
a pre-existing and pre-enriched medium likely by type Ib/Ic 
supernovae (enrichment by type II supernova will involve 
also an O overabundance). 

An account of the absorption can be obtained by study- 
ing the spectral energy distribution (SED) from UV to nIR. 
This is made possible by the detection of GRB 090618 
by Swift/XJVOT and by several robotic telescopes at early 
times. Fitting the SED at 600 s one can infer, in addition 
to the Galactic absorption of Av = 0.29, an intrinsic ab- 
sorption of about Av(z) ~ 0.13 ± 0.06, using a Small Mag- 
ellanic Cloud (SMC) extinction law (x 2 /dof = 6.3/7). A 
fit with a Large Magellanic Cloud or Milky Way extinc- 
tion laws provides unsatisfactory fits (x 2 /dof = 18.9/7 and 
X 2 /dof = 23.2/7). We note that this value is somewhat 
smaller than the value reported by Cano et al. (2011) at 
later times of Av{z) ~ 0.24 ± 0.09, but consistent within 
the errors. 

For a SMC extinction one would expect Nh = Ay x 
1.6 x 10 22 cm 2 (Weingartner & Draine 2001). This im- 
plies for the solar abundance pattern fit with fixed metal- 
licity of the X-ray spectrum a predicted optical absorption 
A v (z) = 0.14 ±0.01 and with free metallicity A v (z) < 0.48. 
Both estimates are in agreement with the observations. If we 



X-ray observations hold a great potential to probe the en- 
vironment surrounding GRBs thanks to the imprint left by 
metals on the afterglow spectrum. Here we exploit an XMM- 
Newton observation of the bright GRB 090618 discovered by 
Swift. The XMM-Newton afterglow spectrum is well fit by 
a power law model absorbed by a (known) Galactic com- 



1 An analysis carried out with free z and metallicity fixed to the 
solar value shows that there are two solutions one at z = 0.5 — 0.7 
and another at z ~ 3.4. This indicates that XMM-Newton data 
alone, in this case, are not sufficient to uniquely identify the GRB 
redshift. 
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Figure 5. Optical and UV Spectral energy distribution data fit. 
Data were taken and shifted according to the overall light curve 
at a common time of 600 s after the burst trigger. The continuous 
line fit refers to an SMC extinction curve; the dashed line fit to 
the LMC and the dot-dashed line to the Milky Way. 



GRBs observed during the prompt event also with XRT, 
because triggered on a precursor event, we have 17 events. 
At the time of the XMM-Newton observation GRB 090618 
is one of the brightest indicating that to improve on the 
results obtained in this letter we would need: i) larger in- 
trinsic column densities (in order to have larger absorption 
signatures). This burst has an intrinsic column density of 
2.3 x 10 21 cm -2 , placing it in the tail of the column den- 
sity distribution (Campana et al. 2010), ii) longer exposure 
times (in order to collect more photons, extending by an- 
other 25 ks the observation would have added about 50% 
more photons), Hi) faster response (in order to collect more 
photons, being on target 1 hr before would increase the to- 
tal number of collected photons by ~ 20%), and iv) lower 
particle background (in order to throw away less exposure 
time and in turn get more photons). Along these lines XMM- 
Newton studies of the ambient medium surrounding GRBs 
is promising. 
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Figure 6. Light curves of all Swift XRT GRBs with a count 
rate larger than 1,000 c s — 1 . GRBs observed with XRT during 
the prompt event were excluded. The unabsorbed X-ray flux is 
in the 0.3-10 keV energy band. The light curve of GRB 090618 
is marked with filled dots (red in the on-line version), all the 
other light curves with dashed lines. The two vertical dashed lines 
indicate the start and the end of the XMM-Newton observation. 



consider instead the fit with a peculiar abundance pattern 
(i.e. the one with Ne, S and Fe free) the intrinsic column 
density is larger and the predicted Av(z) — 0.751qJ is just 
barely consistent with observations (see Fig. 5). We note 
that it is not common that X-ray and optical absorption 
are in agreement. 

Finally, we try to put in a context our exploratory study 
of bright GRBs with XMM-Newton (see also Campana et 
al. 2011). In Fig. 6 we show the light curves of Swift GRBs 
with a count rate larger than 1,000 c s _1 . Excluding those 
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